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Abstract — The paper considers the solute transport in a heterogeneous porous 
medium consisting of well-permeable and poorly permeable zones, taking into 
account nonequilibrium adsorption in the zones. In a well-permeable zone, there are 
two areas, in each of which there is an adsorption of a substance with reversible 
nonequilibrium kinetics. The exchange of solute with the second zone is modeled by 
the source term in the form of a fractional time derivative of the concentration of the 
substance in the first zone. The numerical implementation of the model is carried out 
and the effect of mass transfer to the second zone on the characteristics of the solute 
transport in the first zone is estimated. 


Index Terms: adsorption, adsorption kinetics, approximation, fractional 
derivative, porous medium. 


I. INTRODUCTION 

Underground aquifers and oil and gas strata, as a rule, have a heterogeneous structure at the 
micro and macro scales [1]. Macroscale heterogeneous formations consist of different zones with 
differing, sometimes very strongly, reservoir properties, i.e. porosity, permeability, etc. Areas with 
good porosity and permeability are good conductors for fluids and various substances suspended or 
dissolved in formation fluids. A typical example of heterogeneous formations are fractured-porous 
media (FPM) [2, 3], the structure of which is represented as a system of cracks surrounding porous 
blocks. 

Colloidal particles can move relatively faster and move long distances in structured porous media 
than in media with a homogeneous structure [5,6,7,8]. In such environments, when modeling the 
solute transport, it is usually assumed that the main paths of movement of a fluid and, along with it, 
suspended solids or dissolved substances in it are fractures. Porous blocks in simplified models are 
considered as impermeable for fluid, but particles or solute can diffuse into them. Thus, two zones 
are formed in the medium, one with a mobile fluid (cracks), and the other with a stationary one 
(porous blocks). Mass transfer processes take place between the zones. 

The advance spread of substances in a porous medium can be result of many factors. Therefore, 
there are certain difficulties in mathematical modeling of this phenomenon. Some models in this 
direction were presented in [9,10,11,12]. These models used the two-zone approach noted above. 
Mass transfer between zones is modeled by a first-order kinetic equation [13, 14]. A somewhat 
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different approach, combining kinetic and equilibrium mass transfer between zones, was proposed 
in [15]. Some modification of the two-zone approach is an approach that takes into account the 
movement of fluid in both zones, but with different scales [9, 10]. 

In [20], the process of transfer of colloids was considered under the assumption that reversible 
retention of colloidal particles with different characteristics (parameters) occurs in both zones. 
Relatively speaking, in each zone there are two regions, in one of which the rate of retention of 
particles is higher than in the second, and the rate of release of particles is relatively low than in the 
second. This approach is more realistic than in [4], since irreversible retention of particles is 
observed only at the initial stage of the process, when a monolayer of colloidal particles is formed 
on the surface of the medium rock, i.e. when there is no energy barrier for retention [17,18,19]. 
Thus, the irreversible retention of particles has a time limit, beyond which it becomes reversible. In 
[4], irreversible retention is used for the entire duration of the process, i.e. for the entire time range. 
In this regard, the use of double reversible retention kinetics is preferred. In this case, it becomes 
possible to consider the irreversible kinetics of retention as the limiting case of reversible kinetics. 

In [15], the solute transport in a medium with double porosity is considered, taking into account 
the reversible and irreversible deposition of matter in both zones and the first-order equilibrium 
exchange between the zones. In each zone, i.e. in fractures and porous blocks, there is a reversible 
and irreversible sedimentation (deposition) of a substance with different characteristics, described 
by linear equations. An analytical solution to the problem was obtained, which was used to describe 
the results of earlier experiments [16]. 

In this work, in contrast to [15], a new model is proposed, where the presence of the second zone 
of an inhomogeneous medium is taken into account in the form of a sink (source) term in the 
transport equation written for the first zone. The stock term is presented as a fractional time 
derivative of the concentration of the substance in the first zone with a certain coefficient. This 
model is implemented numerically. The solution is compared with the solution [15] in the first zone. 
It is shown that by appropriate choice of model parameters the solution can be approximated to the 
solution [15]. 


I. MATHEMATICAL MODEL 

A heterogeneous porous medium is considered, consisting of well-permeable and relatively 
poorly permeable zones, the diagram of which is shown in Fig. 1. The parameters in the first zone 
are indicated by the index 1. There are two sections in zone 1, in each of which the precipitation of 
a substance with reversible nonequilibrium linear kinetics occurs. It is believed that such processes 
also occur in the second zone, but we will not write equations and conditions for it. Substance 
exchange occurs with the second zone, which will be modeled by the fractional time derivative of 
the substance concentration in the first zone [20]. Consequently, in contrast to [15], the 
concentration field in the second zone is not considered. 
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Fig. 1. Scheme of material transfer in a two-zone medium 
The equations of solute transport in the one-dimensional case are written in the form 


041 41 Ocy OY cy = 01 Ocy 
P t +p at +6: at + a2 at = 0D, ax2 Ov ax’ (4) 


where a, — is the coefficient due to the solute transport into the second medium, 867 1, y —is the 
order of the derivative. 
The sedimentation of matter in each of the sections of the zones occurs reversibly in accordance 


with the kinetic equations 
OSa1 
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where 64, ks1 — coefficients of deposition of matter from the fluid phase to the solid phase, 5 
01: 501 —the coefficients of separation of the substance from the solid phase and the transition 
to the fluid, s~?. 

Let a fluid with a constant concentration of a substance be pumped into a medium initially 
saturated with a pure (without substance) fluid from the initial moment of time cg. Let us consider 
such periods of time where the concentration field does not reach the right boundary of the medium 
x — oo, Under these assumptions, the initial and boundary conditions for the problem have the form 


C1(0, x) z 0, Sai (0, x) - 0, S51 (0, x) = 0, (7) 
C1(€,0) = Co, (8) 
ðc - 
ae te co) = 0. (9) 


IT NUMERICAL APPROXIMATION MODEL 


Problem (4) — (9) is solved by the finite difference method [21]. In the considered area Q = 
{(t,x),0 <t <T,0 > x > œ} introduced a grid uniform in directions 
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where I —is a sufficiently large integer selected so that segment [0, x], x; = ih, overlaps the area of 
the calculated change in fields 01,501 and S;;,h — the grid step in the direction of x, z — is the 
grid step in time +. 

In the open grid area 


(t; xi); 5 = Tj, Xi = ih, 
Wrh = 7 oo te 
m re ASISTIA] 
equations (4), (5), (6) were approximated as follows 
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where (c4 j, Sai j, Si J _ grid function values c(1,x), Sa1(t, x), Salt, x), at the 
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point (t;, xi). 


From explicit grid equations (11), (12) we determine (Sai)! (Ss1)/ ₪ 
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Grid equations (10) are reduced to the form 
Alei - ו‎ + Be Gna = —(F,)j, (15) 
where 
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The following procedure for calculating the solution is established. By (13), (14) (S ea) 1 


(S D ** are determined, then the system of linear equations (15) is solved by the Thomas’ 


algorithm with respect to Gy Since pp1, qp1 > 1, schemes (13), (14) are stable, and for (15) 
the stability conditions of the Thomas’ algorithm are satisfied. 


IV. NUMERICAL RESULTS AND THEIR ANALYSIS 


Numerical experiments were carried out according to the above algorithm. 
The calculations used Le following values of the initial parameters: 
5 v, = 75 , ,יט = ן2‎ ₪ = 0.005% 0, = 0. a = 6:1074, kai = 3. 
1074 s7t, kadı = 25. 107 8 red 10s ke = 210s p = 1800 kg/m? 
and various y. 

Some typical results are shown in Fig. 2. As seen from Fig. 2, the outflow of the substance into 
the second zone leads to a delayed propagation of the concentration profiles of the substance in the 
mobile fluid. As a consequence of this phenomenon, delays are also observed in the concentrations 
of the adsorbed substance. 
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Fig. 2. Concentration profiles €, (a), Sa ©), Sp (c) at 
t = 3600 c, 


y= 0,5 feee ), 0,7 = = =), 0,9 + ==), 1,0( ). 
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5. CONCLUSION 


In this work, in contrast to [15] a new model is proposed, where the presence of the second zone 
of an inhomogeneous two—zone medium is taken into account in the form of a sink (source) term in 
the transport equation written for the first zone. The stock term is presented as a fractional time 
derivative of the concentration of the substance in the first zone with a certain coefficient. Thus, this 
approach is mono-continuous, while in [15] the bicontinuous approach was used. 

The transport model is analyzed numerically. It is shown that with a decrease in the index from 
unity with the remaining parameters unchanged, the precipitation of matter intensifies. As a result, 
a lag occurs in the development of the distribution of the concentration of a substance in a mobile 
fluid. 
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